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1. Comparing Landau gauge fixing methods 

The formulation of QCD on the lattice does not require gauge fixing. However, for example, 
if one wishes to study the propagators of the fundamental fields, one must choose a given gauge. 
A possible choice is the so called Landau gauge, d^A^ = 0. On the lattice, the Landau gauge is 
implemented by searching for stationary points of 

Fu[g]=Y,Rc{Tr[g{x)U^{x)gHx + fl)]} (1.1) 

over the gauge orbit of U^. As any other local continuous gauge fixing condition ^ §], Landau 
gauge suffers from the so called Gribov copies problem |^], i.e. there are multiple solutions for the 
gauge condition in each gauge orbit. This rises the question on the possibility of a non-perturbative 
definition of the Landau gauge. A solution to this problem is to identify Landau gauge as the search 
of the unique [Q], up to a global gauge transformation, absolute maxima of < \l.l\ ). Unfortunately, 
the search of a global maximum of a function on a multidimensional manifold is far from being 
trivial. For the particular case of the gauge under discussion, some methods have been devised 
which, hopefully, will be able to compute the absolute maximum of Fu [g] . 

In this work we compare two of them, namely the CEASD method, described in and the 
smeared gauge fixing, described in ||^. The CEASD method combines an evolutionary algorithm 
with the steepest descent (SD) method [Q]. On the other hand, smeared gauge fixing relies on the 
smooting of the gauge fixing hypersurface by smearing the configuration. 

We tested these two methods on SU(3) 16^ pure gauge Wilson action configurations ^ for 
j3 € {5.7,5.8,5.9,6.0,6.2}. For each j8 value, five configurations were generated. For each 
configuration, 500 SD gauge fixings were performed, starting from different randomly chosen 
points. This procedure gives an idea on the number of Gribov copies for each configuration, and 



defines a candidate for the global maximum of (1.1). As expected, the number of Gribov copies 



seems to increase with the physical volume of the lattice. 

For each configuration, we have compared the global maximum obtained with the different 
500 SD ('GMAX' in the figures) with the maxima computed by the CEASD method ('CEASD' 
in the figures), by smeared gauge fixing ('SMGF' in the figures) and by a single steepest descent 
applied to the configuration ('SD' in the figures). 

The various plots in figure [T] show that the smeared gauge fixing only once was able to iden- 
tify correctly the global maximum (j3 =5.9 configuration nr. 6000). On the other hand, the 
CEASD method was successful in all j8 = 6.0 and j8 = 6.2 configurations, and in 4 of the j8 = 5.9 
configurations. Curiously, the CEASD method found a larger Fu than the one found by 500 SD for 
one of the j8 = 5.8 configurations. In what concerns the computation of the absolute maximum of 
Fu [g] , the CEASD method seems to be superior, but it is very computationally demanding. Indeed, 
it requires ~ 30 hours in a single Pentium IV, to be compared with 3 hours for the smeared gauge 
fixing method or ~ 30 minutes for the steepest descent. 



'All the configurations used in this work were generated with MILC code |jl^. 
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Figure 1: Fy maxima found by the different gauge fixing methods considered. 
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2. Gribov copies and the gluon and ghost propagators 



The investigation of the Gribov copies effects on the gluon and ghost propagators has been 
an active research field for some time — see for example, [||, ^ 11 ] and references there in. 



Although in a recent work [ |12[ ] it has been argued that the different maxima of do not change 
the infinite volume expectation values, this is not necessarily true on a finite lattice. 



In this section, the effects of Gribov copies in the gluon [13] and ghost propagators are inves- 
tigated using 302 configurations for a 12'* lattice with j8 = 5.8. For each configuration, 500 SD, 
starting from different randomly chosen points, were performed and saved the gauge configurations 
associated to the largest maximum of Fu, the smallest maximum of Fu, and a few more random 
gauge configurations associated to intermediate values of Fu maxima. Then, the propagators com- 
puted have been labelled by (Fu), the mean value of the corresponding Fu associated to each gauge 
fixed configuration. 

In this work, we compute the gluon propagator D{q^) using the following definitions (see [13] 
for more details) 
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D{q^) = ^ :—T {Tr:\AJq)AJ-q) 



and 



where 



and 



D{q^=0) 



f y 



Tr[A^iq)A^i-q) 



^ju(^) =52 exp[-/^(x + a^^/2)] A^(x + a^^/2), 



A^{x + aef^/2) 



1 



U^{x)-UUx) 



1 



-Tr 



U^{x)-Ul{x) 



On the lattice, the discrete momenta available are given by 

271?!,, 



^11 



aL 



0, l,...L^/2, 



and 



2 . fq^ia 
<7^ = - sm — 

a \ 2 

The ghost propagator was computed using a plane wave source 



and the ghost scalar function computed from 

G(^2^ _ 



-yIg-(^). 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



In figure ^ we report the dependence of the gluon propagator in (Fu). The data shows that 
D{0) increases as (Fy) increases. On the other hand, the low non-zero momenta D{q^) seem to be 
decreasing as (Fu) increases. The high momenta region seems to be blind to the choice of copy. 



4 



Gauge fixing methods and Gribov copies effects in lattice QCD 



Paulo J. Silva 



37.5 r 

37 ■ 
36.5 
36 

la. 35.5 
35 - 

34.5 ■ 
34 

33.5 ■ 



0.8386 0.3338 0.8390 0.8392 0.8394 0.8396 0.8398 0.8400 



(a) n^ = (0,0,0,0) 



tt. 13.6 



0.3386 0.8383 0.8390 0.3392 0.8394 0.8396 0.3398 0.8400 



(b) n^ = (1,0,0,0) 




0.8386 0.3338 0.8390 0.8392 0.8394 0.8396 0.8393 0.8400 

(Fu) 



(C) (1,1,0,0) 



4.2 - 



0.3386 0.8383 0.8390 0.3392 0.8394 0.8396 0.3398 0.8400 



(Fu) 

(d) n^ = (l,l,l,0) 



2.85 
2.3 
2.75 

2.7 



0.8386 0.8388 0.8390 0.8392 0.8394 0.8396 0.8398 0.8400 
{Fu) 



0.0925 r 
0.092 
O0915 

0.091 - 

0.0905 ■ 

^ 0.09 - 

q" 0.0895 ■ 
0.089 

0.0885 - 

0.088 ■ 
0.0875 

0.087 ■ 



0.8386 0.8388 0.8390 0.8392 0.8394 0.8396 0.8398 0.8400 
{Fu) 



(e)n^ = (l,l,l,l) 



(f) n^ = (6,6,6,6) 



Figure 2: Gluon propagator veriw* (Fu). 
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Figure 3: Ghost propagator veraM* {Fu). 
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These results show that Gribov copies are important only for the infrared regime. In what 
concerns the ghost propagator, see figure their values decrease as (Fu) increases. Once more, 
the effect is bigger for low momenta (a 20% effect for the lowest momentum, reduced to 0.2% for 
the highest momentum considered). 

3. Conclusions 

For the lattice Landau gauge, we have compared two methods of gauge fixing tunned to ap- 
proach the global maximum of Fu[g]. For a 16'* lattice, the CEASD method seems to perform 
better than smeared gauge fixing. Unfortunately, CEASD is quite computational demanding, not 
only in what concerns memory but also in CPU time. 

In the second part of this paper, we investigated the effect of Gribov copies in the gluon and 
ghost propagators. Our data shows a clear dependence of gluon and ghost propagators on (Fu). 
Furthermore, the data shows that the changes on the propagators due to the choice of the maximum 
of Fu [g] are mainly on the infrared region. 
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